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Abstract 

The statistical methods based on the classical Boltzmann-Gibbs (BG) approach are at 
heart of essentially all descriptions of multiparticle production processes. In many cases, 
however, one observes some deviations from the expected behaviour. It is also known that 
conditions necessary for the BG statistics to apply are usually satisfied only approximately. 
Two attitudes are possible in such situations: either to abandon statistical approach trying 
some other model or to generalise it to the so called nonextensive statistics (widely used in the 
similar circumstances in many other branches of physics). We shall provide here an overview 
of possible imprints of non-extensitivity existing both in high energy cosmic ray physics and 
in multiparticle production processes in hadronic collisions, in particular in heavy ion collisions. 

PACS numbers: 05.40.Fb 24.60.-k 05.10.Gg 
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1 Introduction 



The high energy collisions are usually connected with production of large number of secondaries 
(mostly TT and K mesons). The strong interactions involved here make their detail descriprion 
from first principles impossible and one is forced to turn to phenomenological models of various 
kinds. The statistical models were the first successful approaches to the multiparticle production 
processes since the beginning of the subject almost half century ago jlj and they remain still very 
much alive today, especially in all analysis of multiparticle data performed from the point of view 
of the possible formation of the new state of matter - the Quark Gluon Plasma (QGP) Q. They 
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are all based on the Boltzmann-Gibbs form of entropy, which is identical to the so called Shannon 
entropy used in the information theory approach. 



We would like to stress at this point that information theory can be also applied to hadronic 
processes Its importance there is best visualized by the following example. Suppose that some 
experiment provides us with data on a limited set of n observables: Ri, R2, ■ . ■ , Rn- This triggers 
the theory, which follows with a variety of models, usually completely different (if not contradictory) 
but each claiming good agreement with these data. The models provide therefore apparently 
different views of the same physical reality (data) leading to confusion on what is actually going 
on. It happens this way because experimental data contain only limited amount of information to 
which theoretical models add their own specific assumptions born not in data but in our minds. To 
quantify this situation one must find the way to answer (in a model independent way) the following 
question: given physical assumptions of a model (plus phase space and conservation laws), what 
are their most trivial consequences? This can be done only be by resorting to information theory 
which defines triviality as the lack of information (the less information the more trivial). The 
measure of information provided by Shannon entropy (or some other kind of entropy, in our talk 
this will be Tsallis entropy Q) allows us then to quantify the whole problem. 

Let us elaborate this point a bit further. Experiment gives us usually set of points to which 
one tries to fit some distribution function p{x), which can be thought to be some probability 
distribution. In this case Shannon entropy corresponding to it is 

p{x) ln[p{x)]. (1) 

We are looking for the most trivial p{x), i.e., for one containing the least information, and it means 
for such, which informational entropy is maximal under given circumstances (constraints)^]. What 
are those constraints? At first, being a probability distribution p{x) must be non-negative and 
normalized, i.e., 

p{x)dx = 1. (2) 



Further, p(x) has to reproduce the experimental results given in the form of expectation values 
(^(fc)(2;)^ — of measured set of observables R^''\ i.e., 

dxR'-''\x)p{x) ^ Rk] k = l,...,m. (3) 

Using eqs.(^) and (^ as constraints one then maximalizes information entropy (|l|) and as results 
one gets the seek for distribution p{x), 



L fc=l 



(4) 



where Z is calculated from the normalization condition. This distribution reproduces known in- 
formation (i.e., it tells us the truth, the whole truth about experiment) and contains the least 
information (i.e., it tells us nothing hut the truth). In this sense it is the unique, most plau- 
sible (least biased), model independent probability distribution describing the outcomes of our 

^One should remember that formula (jl|) assumes that p(x) is measured against some other distribution po(x) 
corresponding to equal a priori probability being assigned to equal volumes of phase space. This is sometimes called 
"the first principle of statistical mechanics" f^. The Boltzmann constant, which normally appears in has been 
henceforth assumed unity, kg = 1. 
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experiment. The Lagrange multipliers /3fc are given as solutions of the constraint equations and 
Ru=d\nZ{(3^,...,P^)/dh- 

Now we can understand why apparently completely disparate models can describe sucessfully 
the same experiment. It is because they share (usually in an implicite than explicite way) the 
same set of basic common assumptions (i.e., constraints Rk) and, at the same time, they differ 
completely in the rest of them, i.e., in those Rk, which contain assumptions particular for a given 
model (and which are completely not relevant to its agreement with experiment)^. 

This example tells us that statistical approach is not necessarily connected with any thermo- 
dynamical model but has, in fact, much broader range of applications. In our example above 
the formula ^ is identical with formulas used in thermal models ||] but parameters (3k can take 
any values depending on the energy and multiplicity of the event under considerations |^ . It is 
worth to remind at this point the very old observation |^ that exponential form (^ of so many 
distributions observed experimentally can have its origin in the fact that the actual measurements 
extend only to subsystems of a total system. The summation over the unmeasured (i.e., essentially 
overaged over) degrees of freedom introduces the randomnization, which can be most economically 
described as a kind of thermal bath leading to thermal-like exponential distributions with some 
effective "temperature "T. 

In this approach dynamics lies in our ability to reproduce the known information (i.e., is hidden 
in (3k)- However, the big advantage of this method is that it prevents us from smuggling in anything 
additional, not being present in data. The outcome depends on the funtional form of the entropy. 
The Shannon form (|l|) is the most popular one because it coincides in thermodynamical limit with 
the Boltzmann-Gibbs (BG) form allowing for contact with standard thermodynamics. However, as 
newtonian mechanics is not eternal and universal (in the sense that in the limit of large velocities it 
must be replaced by special relativity, in the limit of vanishingly small masses quantum mechanical 
effects must be accounted for and in the limit of extremaly large masses (stars, galaxies,...) it must 
be replaced by general relativity) the same can be said about the classical BG statistical mechanics. 
It is fully satisfactory only if 

• effective microscopic interactions are short-ranged (i.e., we have close spacial connections); 

• effective microscopic memory is short-ranged (i.e, we have close time connections); 

• boundary conditions are non (mult i) fractal (i.e., the relevant space-time and/or phase space 
is non (multi)fractal). 

As a matter of fact despite the popularity of BG it seems that in majority of situations we are 
dealing with, at least one of these conditions is not satisfied. Therefore it is justified to look for the 
possible generalization of BG which would provide a compactifications of all possible deviations 
from the ideal BG statistics into a minimal - possible only one - additional parameters. In this 

■^A completely forgotten example of such situation is discovery Q that to explain all experimental data on the 
multiparticle reactions of that time it is only necessary to assume that the phase space is effectively one-dimensional 
(i.e., the px cut-off) and that not the whole initial energy yfs of reaction is used for the particle production (i.e., 
the existence of the inelasticity K or of the leading particle effect). 
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way we arrive at the notion of Tsallis statistics based on Tsallis entropy Q : 

5,^-1^1^^. (5) 
1-9 

We shall not dwell too much on it here as this is the subject of separate lecture by Prof. Tsallis 
himself |^]. The only points we would like to stress here for our purposes are: 

• for (7 — !■ 1 one recovers the usual Shannon entropy form of 

• with the normalization condition imposed on p[x) as before and with the new definition of 
the mean values, {R)q — J dxR{x)p{xY , one replaces formula (|^) with 



(6) 



l-(l-g)^/3,i?W(a;) ; 
fc=i 

(notice that, again, for q ^ 1 (^) goes to (^). As before, this distribution reproduces known 
information and conveys least information but this time information is measured in a different 
way, which is paramatrized by the parameter 

• the resultant staistics is non-extensive in the sense that for two subsets A and B such that 
probabilities Pij{A + B) = pi{A)pk{B) the entropy is not additive but instead 

Sg{A + B) = Sg{A) + Sg{B) + (l - q)Sg{A)Sg{B) (7) 

and additivity is recovered only in the limit g — )■ 1; 

• the whole Legendre structure of thermodynamics is preserved for 5 ^ 1 |Q, . 

Although the nonextensive statistics were used for a long time in many branches of physics |^ , 
their possible imprints in multiparticle production are very fresh. In this talk we shall concentrate 
on some cosmic ray experiments and accelerator experiments, which provided us with examples of a 
nonextensive behaviour. Both will demonstrate traces of nonextensivity seen as nonexponentiality 
of some distributions, i.e., where exp(x) exp (q) = [1 - (1 - At the end we shall, 

however, list all other works in which one can find similar conclusions 



2 Cosmic Ray Example 

Part of the cosmic ray detector system consists of extrathick lead emulsion chambers. Such chamber 
is build with a number of sandwiches of X-ray sensitive film, plate of lead and layer of emulsion 
(first allows to localize the point at which external particle of cosmic rays enters, second forces it to 
interact and to start the electromagnetic cascade process, and third allows to calculate deposited 
energy). In some experiments one observes that distribution of the points of first interaction 
(cascade starting points) does not follow the expected simple exponential rulep|pOt| : 

— = const - exp - - (8) 



dT " V A 



^There is a formalism, which expresses both the Tsallis entropy and the expectation value using the escort 
probability distribution, Pi = p\l X/P^Q' However, it is known that such an approach is different from the 
present nonextensive formalism since the Tsallis entropy expressed in terms of the escort probability distribution 
has a difficulty with the property of concavity. 

*Here A'^ is the number of cascades originating at depth T expressed in the so called cascade units where 1 c.u. 
= 6.4 g/cm-^ = 0.56cm of Pb and A is the mean free path for the p-Pb interacions. 
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but deviates from it noticeably for larger T. This phenomenon acquired even its own name: long 
flying component (because original hadrons tend to fly longer without interaction than naively 
expected). 
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Fig. 1 Depth distribution of the starting points, dN(T)/dT, of cascades in Pamir lead chamber. 
Notice the non-exponential behaviour of data points (for their origin cf. which can be 

fitted by Tsallis distribution (10) with q = 1.3. (This figure is reproduced from Fig. 1 of 
0)- 



In |lO[| we have argued that the observed effect can be just another manifestation of the fiuctu- 
ation of the corresponding hadronic cross section a — Ara^ \ (where A denotes the mass number 
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of the target and rriN is the mass of the nucleon, such a possibihty is widely discussed in the 
hterature and observed in diffraction dissociation experiments on accelerators, cf. for details 
and references). As was shown in such fluctuations of cross section (i.e., in effect 1/A) with 
relative variance 

u; = — > 0.2 (9) 

allow to describe the observed effect. However, it turns out jllj that the same data can be fitted 
by the nonextensive formula 



dN 

= const ' 

dT 



T 



(10) 

with parameter q = 1.3 (in both cases A = 18.85 ± 0.66 in c.u. defined above), cf. Fig. 1. 



This example provides us hint on the possible physical meaning of the nonextensvity parameter 
q. Comparing both explanations one can suspect that perhaps, at least in this case, it has something 
to do with the fluctuations of the seemingly constant parameter 1/A in the exponential formula (^). 
We shall argue later that this is indeed the case but before let us demonstrate another example of 
this kind, now from the high energy multiparticle production processes domain. 



3 Heavy ion collision example 

Heavy ion collisions are of particular interest at present because they are the only place were one 
expects the possible formation of the QGP - a new state of matter [Q. The termal and/or all kind 
of statistical models are playing very important role here because one of the crucial parameter 
is the temperature T of reaction deduced usualy from the transverse momentum distributions, 
dN/dpT- This is, however, valid procedure only when dN/dpT is of exponential form. Therefore 
any deviation from such behaviour are always under detailed scrutiny in which one is searching 
for the possible causes. In jl^ it was suggested that the extreme conditions of high density and 
temperature occuring in ultrarelativistic heavy ion collisions can lead to memory effects and long- 
range colour interactions and to the presence of non-Markovian processes in the corresponding 
kinetic equations jl^]. Actually, as has been shown in jl^, dN{pT)/dpT are best described, cf.. 
Fig. 2, by a slightly nonexponential function of the type 



— ; = const ■ 

dpT 



1-(1-.)^^ 



const.exp . (11) 



Here m is the mass of produced particle, k is the Boltzmann constant (which we shall, in what 
follows, put equal unity) and T is, for the q ^ 1 case, the temperature of the reaction considered 
(or, rather, the temperature of the hadronic system produced) . The deviation from the exponential 
form one finds is very small, on the level q = 1.015. However, as we shall demonstrate later it 
can lead to a quite dramatic effects. It was also shown in that to first order in \q — 1| the 
generalized slope becomes the quantity 

Tg + {q- l)mT. (12) 

with T being temperature of a purely thermal source. This should be contrasted with the empirical 
relation for the slope parameter T, from which the freeze-out temperature (at which hadrons are 
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created from the QGP) T/ is then deduced, 

T = Tf + m{v±)^. 



(13) 




Fig. 2 The results for pr distribution dN{pT)/dpT- notice that q = 1.015 results describes also 
the tail of distribution not fitted by the conventional exponent (i.e., q = 1). This figure is 
reproduced from Fig. 3 of . 

The {v±) is a fit parameter usually identified with the average collective (transverse) flow velocity 
of the hadrons being produced. In ( p^ ) one has, instead, a purely thermal source experiencing 
a kind of blue shift at high mx (actually increasing with mx). The nonextensivity parameter q 
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accounts here for all possibilities one can find in and could therefore be regarded as a new 
way of presenting experimental results with g 7^ 1, signaling that there is something going on in 
the collision that prevents it from being exactly thermal-like in the ordinary sense mentioned above. 

To the same cathegory belongs also analysis of the equilibrium distribution of heavy quarks 
in Fokker-Planck dynamics expected in heavy ion collisions Not going into details in this 

case, we shall only mention that it was found that thermalization of charmed quarks in a QGP 
surroundings proceeding via collisions with light quarks and giuons results in a spectral shape, 
which can be described only in q statistics with q=1.114 (for mc = 1.5 GeV and temperature of 
thermal giuons Tg = 500 GeV). 



4 Fluctuations as a cause of nonextensivity? 

Guided by the cosmic ray example we would like to discuss now the hypothesis that the examples 
of nonextensivity in which exponential distribution becomes a power-like (Levy type) distributions 
is caused by the fluctuations in the system which were not accounted for and which make q ^ 1 
pq, M. In other words, q summarizes the action of averaging over 1/A: 



exp 



1 - 

'^0 



(14) 



where I/Aq = (1/A). This connection, which is clear in the cosmic ray example, should also be 
true for the heavy ion collision. Actually, this later case is even more important and interesting 
because of the long and still vivid discussion on the possible dynamics of temperature fluctuations 
and because of its connection with the problem of QGP production in heavy ion collisions p^ . 
Actually the conjecture ( pl| ) is known in the literature as Hilhorst integral formula |19 but was 
never used in the present context. It says that (for q > I case, where e G (0, 00)) 



e 1 
Ao a 



1 



Via) 



drj jf^ ^ exp 



Ao a 



where a 



9-1 



Changing variables under the integral to r\ 



one obtains 



^g>i(£;Ao) — C9 1 + T 

Aq a 



= C„ 



where /(1/A) is given by the following gamma distribution: 



9>1 



(15) 

(16) 
(17) 
(18) 

Notice that, with increasing a the variance (|18| ) decreases and asymptotically (for a 00, i.e, for 
q ^ I) the gamma distribution ( p^ ) becomes a delta function, /q>i(l/A) = (5(A — Aq). The relative 
variance for this distribution is given by 



A'Ao 



r(a) 



(aJ ^"H-a 



with fi — aXo and with mean value and variation in the form: 

^\ = ± and Hi]' 




aXl' 



1 

= - = q-l. 
a 



(19) 
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For the q < 1 case e is limited to e S [0, Ao/(l — q)]- Proceeding in the same way as before, i.e. 
making use of the following representation of the Euler gamma function (where a' — ~a — jh^) 



a'Xo 



a'Ao - 



1 



r(a') 



dr] ^ exp 



a'Ao 



(20) 



and changing variables under the integral to 77 — — j- 
(Hq) but with a. —a' and with the respective /(1/A) 



we obtain Lq^i{e; Aq) in the form of eq. 



/g<i(l/A) given now by the same gamma 
distribution as in (17) but this time with a — > a' and fi — /i(e) = a'Ao — £• Contrary to the q > 1 
case, this time the fluctuations depend on the value of the variable in question, i.e., the mean value 
and variance are now both e-dependent: 



1 



Ao 



and 




1 



a' 



1 



(21) 



However, the relative variance 



(if 



1 



= 1 



(22) 



remains e-independent and depends only on the parameter q. As above the resulting gamma dis- 
tribution becomes a delta function, /q<i(l/A) = S{X — Aq), for a' — > 00, i.e., for q ^ 1. 



Summarizing: one can indeed say that the nonextensivity parameter q in the Lq{e) distributions 
can be expressed by the relative variance uj of fluctuations of the parameter 1/A in the distribution 



lie)-- 



(J = 1 ± w for q> I {+) and g < 1 (-). 



(23) 



Actually the above result (^3|) is derived for the particular (Gamma-like) shape of the fluctuation 
distribution function. How can it be realized? To answer this question let us notice |p^ , ^ that in 
the case when stochastic variable A is decribed by the usual Langevin equation 



dX 
It 



- + m 



X — (f) ~ const > 0, 



(24) 



(with damping constant r and source term (/)g<i — (Aq — ^) ov (f) — (jjq^i = ^), then for the 
stochastic processes defined by the white gaussian noise form of $(i)0 one obtains the following 
Fokker-Plank equation |M for the distribution function of the variable A: 



df{X) 
dt 



with intensity coefficients K12 defined by eq.(E4) and equal to (cf., pTj): 



i^i(A) = (t>- 



X 



DX and K2{X) = 2 D X^ 



(25) 



(26) 



^It means that ensemble mean (^(t)) = and correlator (for sufficiently fast changes) (^(t) S,{t + At)) = 
2D<5(At). Constants r and D define, respectively, the mean time for changes and their variance by means of the 
following conditions: (A(t)) = Aq exp (—7) and (A^(t = 00)) = ^ Dt. Thermodynamical equilibrium is assumed 
here (i.e., t >> r, in which case the influence of the initial condition Aq vanishes and the mean squared of A has 
value corresponding to the state of equilibrium). 
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Its stationary solution has precisely the Gamma-like form we are looking for: 

with the constant c defined by the normalization condition, Jp°° d{l/X)f{l/X) = 1. It depends on 
two parameters: 

Ke) = ^ and = (28) 

with (pq — (/)g>i.5<i and aq — {a, a') for, respectively, q > I and q < 1. Therefore eq. ( p3[ ) with 
oj = is, indeed, a sound possibility with (in the case discussed above) parameter of nonexten- 
sivity q given by the parameter D and by the damping constant t describing the white noise. 



/(A) 



K2iX) 



exp 



.,A'i(A') 



MX') 



5 Back to heavy ion collisions 

With the possibility of such interpretation of parameter q in mind we can now came back to the 
example of heavy ion collisions and see what are its consequences in this case. It is interesting to 
notice that the relatively small departure of q from unity, g — 1 ~ 0.015 |l^, if interpreted in 
terms of the previous section, indicates that rather large relative fluctuations of temperature, of 
the order of AT/T ~ 0.12, exist in nuclear collisions. It could mean therefore that we are dealing 
here with some fluctuations existing in small parts of the system in respect to the whole system 
(according to interpretation of p^) rather than with fluctuations of the event-by-event type in 
which, for large multiplicity N, fluctuations AT/T = 0.06/ VN should be negligibly small fl^ . 
This controversy could be, in principle, settled by detailed analyses of the event-by-event type. 
Already at present energies and nuclear targets (and the more so at the new accelerators for heavy 
ions like RHIC at Brookhaven, now commisioned, and LHC at CERN scheduled to be operational 
in the year 2006) one should be able to check whether the power-like pt distribution dN{pT)/dpT 
occurs already at every event or only after averaging over all events. In the former case we would 
have a clear signal of thermal fluctuations of the type mentioned above. In the latter case one 
would have for each event a fixed T value which would fluctuate from one event to another (most 
probably because different initial conditions are encountered in a given event). 



One point must be clarified, however. The above conjecture rests on the stochastic equation 
(|2^). Can one expect such equation to govern the T fluctuations? To ansewer this question let 
us turn once more to the fluctuations of temperature jlj, ^ ^ discussed before, i.e., to A = T. 
Suppose that we have a thermodynamic system, in a small (mentally separated) part of which the 
temperature fluctuates with AT T. Let X{t) describe stochastic changes of the temperature in 
time. If the mean temperature of the system is (T) — Tq then, as result of fluctuations in some 
small selected region, the actual temperature equals T' — Tq ~ T^{t)T. The inevitable exchange 
of heat between this selected region and the rest of the system leads to the equilibration of the 
temperature and this process is described by the following equation |23[: 



dT 1 
'dt " T 



(T' ~ T) + ng^O, (29) 



which is, indeed, of the type of eq. (24) (here riq<i = and rig>i = 0). This proves the plaus- 
ability of what was said above and makes the event-by-even measurements of such phenomenon 
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very interesting one. 



6 Other imprints of nonextensivity 

The above discussed examples do not exhaust the Ust of the possible imprints of nonextensivity in 
multiparticle production know (or thought of) at present. As it is impossible to review all of them 
here, we shall therefore only mention them. 

In the possible effects of nonextensivity on the mean occupation numbers Uq and on 

the event-by-event fluctuation phenomena have been discussed. It turns out, for example, that 
some charateristics of fluctuations are extremaly sensitive to even small departures of q from unity. 
Such departure can easily mimick the existing correlations or the influence of resonances. 

In an interesting attempt was presented to fit the energy spectra in both the longitudinal 
and transverse momenta of particles produced in the e^e^ annihilation processes at high energies 
using q-statistical model. In this way one can have energy independent temperature T and nonex- 
tensivity parameter q rising quickly with energy from q — 1 to q = 1.2 and reflecting long range 
correlations in the phase space arising in the hadronization process in which quarks and gluons 
combine together forming hadrons (actually, this observation has general validity and applies to 
all production processes discussed here as well). In similar spirit is the work [p5| which attemts to 
generalized the so called Hagedorn model of multiparticle production to g-statistics. 

To the extend to which self-organized criticality (SOC) is connected with nonextensivity 
one should also mention here a very innovative (from the point of view of high energy collision) 
application of the concept of SOC to such processes [p6| . 

The other two examples do not refer to Tsallis thermostatistics directly, nevertheless they are 
connected to it. First is the attempt to study, by using the formalism of quantum groups, the 
Bose-Einstein correlations between identical particles observed in multiparticle reactions p7| , sec- 
ond are works on intermittency and multiparticle distributions using the so called Levy stable 
distributions They belong, in some sense, to the domain of nonextensivity because, as was 

shown in [ p9| , there is close correspondence between the deformation parameter of quantum groups 
used in and the nonextensivity parameter q of Tsallis statistics and there is also connection 
between Tsallis statistics and Levy stable distributions . Some traces of the possible nonexten- 
sive evolution of cascade type hadronization processes were also searched for in | |30| . The quantum 
group approach |27[ |9| could probably be a useful tool when studying delicate problem of in- 
terplay between QGP and hadrons produced from it. It is plausible that description in terms of 
g-deformed bosons (or the use of some kind of interpolating statistics) would lead to more general 
results than the simple use of nonextensive mean occupation numbers < ?i >g discussed above (for 
which the only known practical description is limited to small deviations from nonextensivity only). 

One should mention also attempts to use Levy-type distributions to fit the development of the 



cosmic ray cascades to learn of how many descendat they contain |31 or the explanation of the 



Feynman scaling violation observed in multiparticle distributions in terms of the q parameter 
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7 Summary 



To summarize, evidence is growing in favour of the view that the standart statistical model can be 
enlarged towards the nonstandart statistics and that by including one new parameter q it allows to 
repoduce much broader set of data than it was done so far. It was demonstrated that q is probably 
connected with the intrinsic fluctuations existing in the system under consideration which were 
previously not considered at all^. Let us close with the remarks that, as was shown in [|^, one can 
also try to use power-like distributions discussed here to the new formulation of the quantum field 
theory (for example, in terms of lorenzian rather than gaussian path integrals, accounting in this 
manner for some intrinsic long range correlations impossible to dispose of and extremaly difficult 
to account for by other methods). 

Acknowledgements: One of us (G.W.) would like to thank the organizers of the RANP2000 for 
support making his participation in this workshop possible. Authors are grateful to Prof. S.Abe 
for his valuable comments. 
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Discussion 



Question by K.Werner: Concerning your example dN/dpT'- this seems to he quite under con- 
troljlarge pt is perturbative QCD (power-like) , small pt is soft, so you have a superposition of 
these two processes. There seems to be no need of these additional fluctuations. 
Answer: Your proposition is just one of many around us. You are right that adding two known 
(in their domain of applicability) mechanisms will also fit data. But what we are saying here is 
that, perhaps, there is mechanism which would show itself in the intermediate region and which 
can be described by the nonextensive parameter q. And it can be experimentally checked in event- 
by-event type of analysis, as was presented above. This would answer the question is there a need 
for something or not. Besides, it should be stressed that using only one new parameter q one can 
fit quite large (in [Q the whole) region oipr- 

Question by F.Grassi: Did you try to fit pt distributions for other type of particles than charged? 
Answer: No, at least not yet. 

Question by T.Kodama: If q represents the measure of fluctuations of ensemble, then it would also 
reflect in multiplicity distributions. However, value of q is different for the pr distribution and 
multiplicity distribution. How do you interpret this? 

Answer: It is difficult for me to comment because so far I have not seen such analysis of the multi- 
plicity data. But if things are really as you say then, assuming that everything was done correctly, 
I would argue that multiplicity is global characteristic whereas px is a local one. Therefore they 
are sensitive to different aspects of the underlying dynamics and resulting parameters q could a 
priori be different. One can also say that for pT distribution the situation is more clear as we are 
simple replacing here exp(...) by exp^ (...). Multiphcity distributions depend on q in indirect way 
only. 
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